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The dimensionless parameter = M 2 / (\6n 2 F 2 ), where F is the pion decay constant in the 
chiral limit and M is the pion mass at leading order in the quark mass, is expected to control the 
convergence of chiral perturbation theory applicable to QCD. Here we demonstrate that a strongly 
coupled lattice gauge theory model with the same symmetries as two-flavor QCD but with a much 
lighter cr -resonance is different. Our model allows us to study efficiently the convergence of chiral 
perturbation theory as a function of S, . We first confirm that the leading low energy constants 
appearing in the chiral Lagrangian are the same when calculated from the e-regime and the p- 
regime. However, < 0.002 is necessary before 1-loop chiral perturbation theory predicts the 
data within 1%. However, for > 0.0035 the data begin to deviate qualitatively from 1-loop 
chiral perturbation theory predictions. We argue that this qualitative change is due to the presence 
of a light c-resonance in our model. Our findings may be useful for lattice QCD studies. 
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1. Introduction 

Chiral perturbation theory is an effective field theory that describes the low momentum physics 
of the pions. According to this theory physical quantities can be expressed as a power series in 
the dimensionless parameter = M 2 /(\6n 2 F 2 ) where M 2 = mL/F 2 is the square of the pion 
mass to the leading order in the quark mass m, F is the pion decay constant and £ is the chiral 
condensate both evaluated in the chiral limit. The unknown coefficients of the power series encode 
the dynamics of QCD and are referred to as the low energy constants. One of the important topics 
of research today is to match lattice QCD data to the chiral expansion and compute these low energy 
constants from first principles [jlj ^]. In order to match data with the chiral expansion reliably it is 
important to find the range in % where the predictions will be valid [Q, |J]. The region of validity is 
governed by the properties of resonances in the theory. Can the a -resonance affect the convergence 
of the chiral expansion? 

The a-resonance arises in % — % scattering in a channel with vacuum quantum numbers. In 
the physical world sigma is a broad resonance and is usually ignored in the analysis. Recently, it 
was estimated that M a ~ 440MeV and T a ~ 544MeV [Q]. On the other hand, in lattice QCD, as 
the pion masses increase, the properties of the resonance will clearly change. Indeed recent studies 
do find a strong dependence on the quark mass [§]. It is interesting to ask if this dependence can 
affect the chiral expansion. This question is non-perturbative and difficult to answer within lattice 
QCD currently. On the other hand it may be answerable in models with the same symmetries as 
QCD and such studies may help shed some light on the subject. 

Here we study a QCD-like lattice field theory model which has the same symmetries as two- 
flavor QCD. Our model also contains a parameter which we tune so that it contains a light sigma 
resonance in addition to pions. We then find evidence that indeed chiral perturbation theory breaks 
down when M„ > M a /3. For more details we refer to the published versions of this work in 

2. Model and Observables 

Our model involves two flavors of staggered fermions interacting strongly with abelian gauge 
fields. The action of the model is given by 

where x denotes a lattice site on a 4 + 1 dimensional hyper-cubic lattice L t xL 4 . Here L 4 is the usual 
Euclidean space-time box while L t represents a fictitious temperature direction whose role will be 
discussed below. The two component Grassmann fields, \j7 x and \\f x , represent the two quark (u,d) 
flavors of mass m, and 0^ is the compact U(l) gauge field through which the quarks interact. 
Here /I = 1,2,.., 5 runs over the 4 + 1 directions. The pL = 1 direction will denote the fictitious 
temperature direction, while the remaining directions represent Euclidean space-time. The usual 
staggered fermion phase factors Tfy Y obey the relations: r\l=T and r\f = 1 for i = 2, 3,4, 5. The 
parameter T controls the fictitious temperature. The four fermion coupling c sets the strength of 
the anomaly. As explained in [Q], the above model has the same symmetries as Nf = 2 QCD. 



(2.1) 
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We have developed an efficient cluster algorithm to solve this model and have studied it in the 
£-regime [Q] and the /^-regime [j|]. In our work we fix L t = 2 and c = 0.3. For these parameters the 
temperature T can be tuned so that the model is either in a spontaneously broken phase for T < T c 
or in the symmetric phase for T > T c , where T c = 1.73779(4) 0]. Since the phase transition is 
second order, close to T c the pion decay constant in the chiral limit F is small in lattice units. This 
reduces the lattice artifacts in our model. Further, tuning T close to T c in the low temperature phase 
also guarantees the existence of a light a-resonance. For this reason, we choose to fix T = 1.7 in 
this work. 

We focus on three observables: The vector current susceptibility Y v , the chiral current suscep- 
tibility Y c , and the chiral condensate susceptibility % a . These are defined as 

*v = ^7 ( I (ew) ) ' x ° = h^ {WxWx WyVy) (2 ' 2) 

where TL (x) and 1° (x) denote one of the components of the vector and the chiral current respec- 
tively. For a detailed discussion of our algorithm and observables, we refer the reader to [j^J. 



3. Results 



We first set m = and study the finite size scaling of our observables. This is in the £ -regime 
of chiral perturbation theory. At 1-loop the theory predicts that [^,[l(], |iT, 12] 



Y v = — 1 + 



0.14046 

~m 2 



+ 



[FLY... 



I?L A ( 0.42138 
I , Xa = —r- 1 + + 



(FL) 2 (FL) 4 ... 



(3.1) 



Our data and fits are shown in figure [jj Using the fits we extract the low energy constants F and £ 
and find that F = 0.2327(1), a = 1.91(9) with a % 2 /DOF = 1.2 and £ = 0.4346(2), b = 1.72(11) 
with a x 2 /DOF = 0.2. 

Next we vary the quark mass in the interval 0.0002 <m< 0.01 for lattices in the range 12 < 
L < 32 and thus explore the p-regime of the chiral expansion. Here the 1-loop predictions for Y c , 
Y v and % a are given by [[K], [I2p. 



Y c = (F„) 



(Fx) 



\-2g~ l (LM JZ )^ + ^ : 



2L 



oL 



Xa = ((qq)) 2 L 4 



\-3g~ l (LM JZ )^+^ : 



(3.2a) 



(3.2b) 



(3.2c) 



where M K is the pion mass, F % is the pion decay constant and (qq) is the chiral condensate at a 
given quark mass m. The function g~\ arises due to pions constrained to be inside a periodic box 
and is given by 

£iW= L t^i(^) (3-3) 

11 ,«2 ,"3,147^0 v 
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L L 



Figure 1: Finite size scaling of observables in the e-regime. The left figure shows the current susceptibility 
and the right figure shows the chiral condensate susceptibility both at m = 0. The insets show a larger range 
of data. The solid line is the fit to the eq.(3. 1) while the dashed line is obtained when we set a = b = 0. 

where K\ is a Bessel function of the second kind and n = n \ + n\ + n\ + n\. Fitting our data with 
these predictions we can determine M n ,F n and (qq) as a function of the quark mass m. 

Our data fit well to the above predictions for 0.0002 < m < 0.0035. For m > 0.0035 the 
fits converge only if we exclude almost all the curvature in Y c and % a . In particular, we are not 
sensitive to the g\(X) function for these two observables and the data fit well even to a constant. 
On the other hand Y v continues to fit well for the entire range of data and this can be used to 
extract M n accurately as a one parameter fit. This may be a useful observation even for lattice 
QCD calculations. Thus, we were able to extract F % , M n and (qq) as functions of the quark mass. 
Some of our results are summarized in Table [!} As an illustration, we also show our data at 
m = 0.0065, m = 0.0035, m = 0.001 and m = 0.0002 along with the fits in Fig. |. 

The quark mass dependence of F n , (qq) and M n have also been computed up to 1-loop in 



m 


{m) 




M n 


0.0002 


0.4392(2) 


0.2348(1) 


0.0400(2) 


0.0005 


0.4441(2) 


0.2377(1) 


0.0627(2) 


0.0010 


0.4528(2) 


0.2423(1) 


0.0878(1) 


0.0020 


0.4678(2) 


0.2501(1) 


0.1220(2) 


0.0035 


0.4867(2) 


0.2606(1) 


0.1584(2) 


0.0050 


0.5024(3) 


0.2690(2) 


0.1860(3) 


0.0065 


0.5170(3) 


0.2764(2) 


0.2083(4) 


0.0100 


0.5433(2) 


0.2912(2) 


0.2521(5) 



Table 1: Results for M n , F K and (qq) from fitting Y V ,Y C , and % a as a function of L to the finite-size 1-loop 
chiral perturbation theory. 
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Figure 2: Finite size scaling of Y V ,Y C and % c at m = 0.0002 (top left), m = 0.001 (top right), m = 0.0035 
(bottom left) and m = 0.0065 (bottom right). The solid lines are fits of the data to the expected finite size 
scaling form from chiral perturbation theory while dashed lines are fits to a constant. 



F n = F 



l-£'(log£'-2c F ) 



; (qq) = L 



i--^(i og r-2c E 



• M 2 ■■ 



\ + -^(\og^-2c M ) 



(3.4) 

where cp,cz and cm are higher order low energy constants and are usually defined in the literature 
as c; = log(A,/47z;F). We have performed a combined fit of F K ,(qq) and M n , quoted in Table [T], in 
the region 0.0002 <m< 0.001 to the above three relations. The result is tabulated in Table |2| The 
values of F and Z obtained in the ^-regime agree nicely with the those computed in the £-regime. 

In order to isolate the region where 1-loop corrections are a good description of the data we 
define the following rescaled and subtracted quantities: Rp = F K /F — 1 + ^'log^', = (qq)/L — 



L 


F 




CF 


cm 


X 2 


0.4354(3) 


0.2329(2) 


11.9(3) 


19.3(5) 


39(3) 


1.1 



Table 2: Results from a combined fit of the data in Table 1 to Eqs.(3.6) 
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Figure 3: Rescaled and subtracted quantities defined in the text which go to zero linearly in the region 
where 1-loop chiral perturbation theory is valid. The solid lines are plots of the fits discussed in the text. 
The dashed lines show the linear region for larger values of The "knee" is estimated roughly as the point 
where the two lines cross. 

1 +3§'log£72, and R M = M 2 n /M 2 - 1 - §'log(§')/2 using F = 0.2329 and I = 0.4354 obtained 
from our fits. By definition, the /?'s must be linear in in the region where 1-loop results are valid. 
In Fig. H we plot the R's as a function of t,' . Assuming errors of 1% or less can be tolerated, Fig. || 
shows that 1-loop chiral perturbation theory describes the data for < 0.002. Interestingly, there 
is also an approximately linear region for > 0.006 but with a completely different slope. This 
is shown as the dashed line in Fig. ||. This behavior suggests that chiral perturbation theory begins 
to break down roughly around ps 0.0035, which is the location of the "knee" that separates the 
low and high t,' regions. We will argue below that the a-resonance is responsible for this break 
down. 



4. Discussion and Conclusions 

It has been argued in the context of the 0(4) linear sigma model, that the physics in the sigma 
channel is directly related to the coefficients cp,cz and cm- Perturbative calculations show that 



[13, 14, 1 



7 $7l 2 7 § 

cz = \og(M R /4nF)-- + — , c M = \og(M R /4nF) - - + (4.1) 

6 igR 3 g R 

where M 2 a = M|[l + g R (?>ny/3> - 13)/(167T 2 )]. Here M a is that mass of the a particle and g R is 
the corresponding renormalized coupling, g R = M\/2F 2 . Since our model is close to the 0(4) 
phase transition the linear sigma model should be a reliable description of the physics in the sigma 
channel. Indeed, using cz = 12 we find that M a /F ~ 2 while using cm = 39 we again find that 
Mo/F ~ 2. The fact that these two agree with each other is a confirmation of our belief. Assuming 
Mo/F ~ 2 and setting the scale of our lattice with F = 90MeV we estimate M a ~ 180MeV in our 
model. At ~ 0.0035 we find that M K ~ 60MeV. Hence, we conclude that when M K > M a /3 
chiral perturbation theory begins to break down and the physics is better described by the linear 
sigma model. 
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There are many important differences between our model and QCD. However the most im- 
portant difference is that our model contains a light and perhaps a narrow a -resonance while in 
QCD the a-resonance is expected to be heavier and broader. This is the reason why our low energy 
constants turned out to be much larger than QCD and the convergence of chiral perturbation theory 
was affected. Despite the differences, it is indeed encouraging that our model supports the most 
important expectations of chiral perturbation theory namely, chiral perturbation theory is applica- 
ble in a region of small quark masses and that the properties of resonances play an important role in 
determining this region. In particular, we have found evidence that the properties of a-resonance 
are encoded in the low energy constants that control the chiral logarithms and hence can play an 
important role in determining the region where 1-loop chiral perturbation theory is valid. The res- 
onance properties of course change with the quark masses. Thus, it is safe to assume that chiral 
perturbation theory can only become reliable in the region of the quark mass where the properties 
of the a-resonance (and other resonances) varies little. A rough estimate based on ^ suggests 
that M n < 250MeV may be necessary. Thus, it should not be very surprising that 1-loop chiral 
perturbation theory may be applicable at a few percent accuracy only at realistic pion masses. 
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